Attorney Docket No.: 1 5685.P049 Patent Application 

IN THE UNITED STATES PATENT AND TRADEMARK OFFICE 



TITLE OF THE INVENTION 



Training Sequence with a Random Delay for a Radio Communications 

System 



Inventor 

MlTHAT C. DOGAN 



Prepared by 

Blakely, Sokoloff, Taylor & Zafman llp 
12400 Wilshire Boulevard 

Seventh Floor 
Los Angeles, CA 90025-1026 
(303)740-1980 



express mail certificate of mailing 



"Express Mail" mailing label number: EL609099423US 

Date of Deposit: ■ Sf^ Arrf^Q 

I hereby certify that I am causing this paper or fee to be deposited with the United States Postal Service "Express Mail 
Post Office to Addressee" service on the date indicated above and that this paper or fee has been addressed to the 
Commissioner of Patents and Trademarks, Washington, D. C. 2023 1 

April M. Worlev 

(Typed or printed name of person mailing paper or fee) 

■VJl 

(Signature of person mailing paper or fee) 

^36/ 




(Date signed) 




Training Sequence with a Random Delay for a Radio Communications 

System 

BACKGROUND OF THE INVENTION 

5 Field of the Invention 

The present invention applies to burst timing for radio communications systems 
and, in particular, to selecting from among different amounts of delay for use in sending a 
burst from one station to another. 

10 Description of the Prior Art 

Mobile radio communications systems such as cellular voice and data radio 
systems typically have several base stations in different locations available for use by 
mobile or fixed user terminals, such as cellular telephones or wireless web devices. Each 

15 base station typically is assigned a set of frequencies or channels to use for 

communications with the user terminals. The channels are different from those of 
neighboring base stations in order to avoid interference between neighboring base 
stations. As a result, the user terminals can easily distinguish the transmissions received 
from one base station from the signals received from another. In addition, each base 

20 station can act independently in allocating and using the channel resources assigned to it. 

Such radio communications systems typically include a broadcast channel (BCH). 
The BCH is broadcast to all user terminals whether they are registered on the network or 
not and informs the user terminals about the network. In order to access the network, a 
user terminal normally tunes to and listens to the BCH before accessing the network. It 

25 will then use the information in the BCH to request access to the network. Such a request 
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typically results in an exchange of information about the network using separate control 
and access channels and ends in the user terminal receiving an assignment to a particular 
base station. 

Before the user terminal has received an assignment, its request for access is 
5 typically sent in a timeslot that is shared with other user terminals. This allows for the 
possibility that its request will collide with other requests. If the request is received at the 
same time as the colliding request, it can be difficult to resolve the colliding requests at 
the base station. However, accurate delay timing advance or delay assignments cannot be 
made by the base station until the range to each of the user terminals is known. In 
10 addition, the BCH does not permit assignments to be made as it is the same for all user 
terminals that receive it. 
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BRIEF SUMMARY OF THE INVENTION 
A method and apparatus are provided that allows a delay to be chosen to be 
applied to a user terminal. According to one aspect of the present invention, the 
invention includes receiving a downlink broadcast burst from a base station, determining 
timing for an uplink burst from the received broadcast burst, selecting an amount of delay 
for the uplink burst, and transmitting an uplink burst to the base station using the 
determined timing and the selected amount of delay. 
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BRIEF DESCRIPTION OF THE DRAWINGS 

The present invention is illustrated by way of example, and not by way of 
limitation, in the figures of the accompanying drawings in which like reference numerals 
refer to similar elements and in which: 

Figure 1 is a simplified block diagram of a base station on which an embodiment 
of the invention can be implemented; 

Figure 2 is a block diagram of a remote terminal on which an embodiment of the 
invention can be implemented; 

Figure 3 is a diagram of an exemplary broadcast channel BCH burst; 

Figure 4 is a diagram of an exemplary configuration request CR burst; 

Figure 5 is a diagram of an exemplary training sequence in accordance with the 
present invention; 

Figure 6 is a flow chart of an implementation of aspects of the present invention; 
Figure 7 is a graph of a search for a least squares error relative timing; and 
Figure 8 is a graph of a follow-up search for a least squares error relative timing. 




DETAILED DESCRIPTION OF THE INVENTION 



Base Station Structure 



5 The present invention relates to wireless communication systems and may be a 

fixed-access or mobile-access wireless network using spatial division multiple access 
(SDMA) technology in combination with multiple access systems, such as time division 
multiple access (TDMA), frequency division multiple access (FDMA) and code division 
multiple access (CDMA). Multiple access can be combined with frequency division 
if 10 duplexing (FDD) or time division duplexing (TDD). Figure 1 shows an example of a 

n 1 base station of a wireless communications system or network suitable for implementing 

\J 

fy the present invention. The system or network includes a number of subscriber stations, 

f™ 

fU also referred to as remote terminals or user terminals, such as that shown in Figure 2. 

The base station may be connected to a wide area network (WAN) through its host DSP 
2? 15 23 1 for providing any required data services and connections external to the immediate 
wireless system. To support spatial diversity, a plurality of antennas 103 is used, for 
example four antennas, although other numbers of antennas may be selected. 

^S^et of spatial multiplexing weights for each subscriber station are applied to the 
respective modbl^ted signals to produce spatially multiplexed signals to be transmitted by 
20 the bank of four anterma^. The host DSP 231 produces and maintains spatial signatures 
for each subscriber station forb^ch conventional channel and calculates spatial 
multiplexing and demultiplexing wei^hte using received signal measurements. In this 
manner, the signals from the current active sb^scriber stations, some of which may be 
active on the same conventional channel, are separated and interference and noise 
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suppressecKWhen communicating from the base station to the subscriber stations, an 
optimized mulficlobe antenna radiation pattern tailored to the current active subscriber 
station connections^and interference situation is created. The same spatial signature is 
applied to signals received from the user terminal by the base station in order to resolve 
5 simultaneous signals on thessame channel. Suitable smart antenna technologies for 

achieving such a spatially directed beam are described, for example, in U.S. Patents Nos. 
5,828,658, issued Oct. 27, 1998 to Stttersten et al. and 5,642,353, issued June 24, 1997 to 
Roy, III et al. \ 

The outputs of the antennas are connected to a duplexer switch 107, which in this 

10 TDD system is a time switch. Two possible implementations of switch 107 are as a 

frequency duplexer in a frequency division duplex (FDD) system, and as a time switch in 
a time division duplex (TDD) system. When receiving, the antenna outputs are 
connected via switch 107 to a receiver 205, and are mixed down in analog by RF receiver 
("RX") modules 205 from the carrier frequency to an FM intermediate frequency ("IF"). 

15 This signal then is digitized (sampled) by analog to digital converters ("ADCs") 209. 

Only the real part of the signal is sampled. Final down-converting to baseband is carried 
out digitally. Digital filters can be used to implement the down-converting and the digital 
filtering, the latter using finite impulse response (FIR) filtering techniques. This is shown 
as block 213. The invention can be adapted to suit a wide variety of RF and IF carrier 

20 frequencies and bands. 

There are, in the present example, four down-converted outputs from each 
antenna's digital filter device 213, one per receive timeslot. The particular number of 
timeslots can be varied to suit network needs. While the present example uses four 
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uplink and four downlink timeslots for each TDD frame, desirable results have also been 
achieved with three timeslots for the uplink and downlink in each frame. For each of the 
four receive timeslots, the four down-converted outputs from the four antennas are fed to 
a digital signal processor (DSP) device 217 (hereinafter "timeslot processor") for further 
processing, including calibration, according to one aspect of this invention. Four 
Motorola DSP56303 DSPs can be used as timeslot processors, one per receive timeslot. 
The timeslot processors 217 monitor the received signal power and estimate the 
frequency offset and time alignment. They also determine smart antenna weights for 
each antenna element. These are used in the spatial diversity multiple access scheme to 
determine a signal from a particular remote user and to demodulate the determined signal. 

The output of the timeslot processors 217 is demodulated burst data for each of 
the four receive timeslots. This data is sent to the host DSP processor 231 whose main 
function is to control all elements of the system and interface with the higher level 
processing, which is the processing which deals with what signals are required for 
communications in all the different control and service communication channels defined 
in the system's communication protocol. The host DSP 23 1 can be a Motorola 
DSP56303. In addition, timeslot processors send the determined receive weights for each 
user terminal to the host DSP 23 1 . The host DSP 23 1 maintains state and timing 
information, receives uplink burst data from the timeslot processors 217, and programs 
the timeslot processors 217. In addition it decrypts, descrambles, checks error correcting 
code, and deconstructs bursts of the uplink signals, then formats the uplink signals to be 
sent for higher level processing in other parts of the base station. With respect to the 
other parts of the base station it formats service data and traffic data for further higher 
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processing in the base station, receives downlink messages and traffic data from the other 
parts of the base station, processes the downlink bursts and formats and sends the 
downlink bursts to a transmit controller/modulator, shown as 237. The host DSP also 
manages programming of other components of the base station including the transmit 
5 controller/modulator 237 and the RF timing controller shown as 233. 

The RF timing controller 233 interfaces with the RF system, shown as block 245 
and also produces a number of timing signals that are used by both the RF system and the 
modem. The RF controller 233 reads and transmits power monitoring and control values, 
controls the duplexer 107 and receives timing parameters and other settings for each burst 

10 from the host DSP 231. 

The transmit controller/modulator 237, receives transmit data from the host DSP 
23 1 , four symbols at a time. The transmit controller uses this data to produce analog IF 
outputs which are sent to the RF transmitter (TX) modules 245. Specifically, the 
received data bits are converted into a complex modulated signal, up-converted to an IF 

1 5 frequency, 4-times over-sampled, multiplied by transmit weights obtained from host DSP 
231, and converted via digital to analog converters ("DACs") which are part of transmit 
controller/modulator 237 to analog transmit waveforms. The analog waveforms are sent 
to the transmit modules 245. 

The transmit modules 245 up-convert the signals to the transmission frequency 

20 and amplify the signals. The amplified transmission signal outputs are sent to antennas 
103 via the duplexer/time switch 107. 

User Terminal Structure 
Docket No. 01 5685.P049 
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Figure 2 depicts an example component arrangement in a remote terminal that 
provides data or voice communication. The remote terminal's antenna 45 is connected to 
a duplexer 46 to permit antenna 45 to be used for both transmission and reception. The 
5 antenna can be omni-directional or directional. For optimal performance, the antenna can 
be made up of multiple elements and employ spatial processing as discussed above for 
the base station. In an alternate embodiment, separate receive and transmit antennas are 
used eliminating the need for the duplexer 46. In another alternate embodiment, where 
time division diversity is used, a transmit/receive (TR) switch can be used instead of a 
10 duplexer as is well-known in the art. The duplexer output 47 serves as input to a receiver 
48. The receiver 48 produces a down-converted signal 49 which is the input to a 
demodulator 51. A demodulated received sound or voice signal 67 is input to a speaker 
66. 

The remote terminal has a corresponding transmit chain in which data or voice to 
15 be transmitted is modulated in a modulator 57. The modulated signal to be transmitted 
59, output by the modulator 57, is up-converted and amplified by a transmitter 60, 
producing a transmitter output signal 61 . The transmitter output 61 is then input to the 
duplexer 46 for transmission by the antenna 45. 

The demodulated received data 52 is supplied to a remote terminal central 
20 processing unit 68 (CPU) as is received data before demodulation 50. The remote 
terminal CPU 68 can be implemented with a standard DSP (digital signal processor) 
device such as a Motorola series 56300 DSP. This DSP can also perform the functions of 
the demodulator 51 and the modulator 57. The remote terminal CPU 68 controls the 
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receiver through line 63, the transmitter through line 62, the demodulator through line 52 
and the modulator through line 58. It also communicates with a keyboard 53 through line 
54 and a display 56 through line 55. A microphone 64 and speaker 66 are connected 
through the modulator 57 and the demodulator 5 1 through lines 65 and 66, respectively 
for a voice communications remote terminal. In another embodiment, the microphone 
and speaker are also in direct communication with the CPU to provide voice or data 
communications. 

The remote terminal's voice signal to be transmitted 65 from the microphone 64 is 
input to a modulator 57. Traffic and control data to be transmitted 58 is supplied by the 
remote terminars CPU 68. Control data 58 is transmitted to base stations during 
registration, session initiation and termination as well as during the session as described 
in greater detail below. 

In an alternate embodiment, the speaker 66, and the microphone 64 are replaced 
or augmented by digital interfaces well-known in the art that allow data to be transmitted 
to and from an external data processing device (for example, a computer). In one 
embodiment, the remote terminars CPU is coupled to a standard digital interface such as 
a PCMCIA interface to an external computer and the display, keyboard, microphone and 
speaker are a part of the external computer. The remote terminal's CPU 68 
communicates with these components through the digital interface and the external 
computer's controller. For data only communications, the microphone and speaker can 
be deleted. For voice only communications, the keyboard and display can be deleted. 
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Broadcast Channel (BCH) 

In one embodiment, the system of the present invention is initiated for each user 
terminal or remote terminal from the broadcast channel BCH which is transmitted as a 
5 burst from the base station to all potential user terminals. The BCH burst, unlike the 
traffic channel bursts, is transmitted in all directions where user terminals may be, 
typically omnidirectionally but the specific beam pattern will depend on the network. An 
example of a broadcast burst structure is shown in Figure 3. The BCH communicates 
enough basic information to enable a subsequent exchange of configuration request CR 
10 and configuration message CM between the base station and the user terminal. The BCH 
also provides general frequency offset and timing update information to all user 
terminals. 

Table 1 , below summarizes the content of an example of a BCH burst, as shown 
in Figure 3. 

15 
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Duration 


Contents 


10 jisec 


ramp - up 


272 usee 


frequency correction training symbols f|, f2 ? ..., fi36 


256 |isec 


timing correction training symbols t !? t 2 ,...ti28 


16 jasec 


broadcast preamble n, r 2 ,...r8 


512 (isec 


information symbols h'i, h^.-h^ 


10 (isec 


ramp - down 


14 jasec 


inter-burst guard time 



Table 1 



The frequency and timing correction training symbols can be set according to any 
one of many approaches well-known in the art. They can also be combined, exchanged 
with a synchronization sequence or eliminated. 

The broadcast information symbols are constructed from a 15-bit broadcast 
message which is modulated and coded into a 256 bit sequence. The number of symbols 
as well as the structure and sequence of transmitted bits can be varied to suit a wide 
variety of applications. The broadcast channel information symbols provide the 
information needed for a user terminal to request a configuration message from the base 
station. 

Each broadcast message is mapped into a broadcast burst with the information 
shown below in Table 2. 
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Broadcast Message 


Field 


# of Bits 


BStxPwr 


5 


BSCC 


7 


BSload 


3 


Total 


15 



Table 2 



BStxPwr is the effective isotropic radiated power of the broadcast message. This 
number indicates the power transmitted by the base station taking into account the 
number of amplifiers and diversity antennas available at the base station. 

BSCC is the base station color code, used by the user terminal to select training 
data for uplink bursts and to distinguish broadcasts of different base stations. 

BSload is an indication of the amount of unused capacity the base station has. 

In one embodiment, the network is designed to take maximal advantage of spatial 
division multiple access technologies and particularly smart antenna array signal 
processing. To help maintain reliable spatial channels in an extremely dense frequency 
reuse pattern, the network uses time division duplex TDMA where uplink and downlink 
transmissions are always on the same frequency. In addition, because many user 
terminals are single antenna and transmit and receive omnidirectionally, except for the 
BCH, an uplink burst is always received before a downlink burst needs to be sent. This 
allows downlink bursts to be more accurately spatially directed. An uplink training 
sequence is embedded in every uplink burst to allow for moderately fast frequency 
hopping despite any decorrelation of the spatial channel with frequency. 
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The frequency hopping sequence may be any one of many different sequences 
well-known in the art. In one embodiment, the parameters of the frequency hopping 
scheme are initially unknown to the user terminal. This maximizes the flexibility of the 
network and increases the flexibility of the user terminal. As explained below, the 
5 frequency hopping parameters are transmitted to the user in the CM burst. 

In one embodiment, the BCH channel is shared by all base stations in the wireless 
communication system. Using the 7 bit BSCC, up to 128 base stations can be 
accommodated. The BCH is part of a time division duplex channel with a repeating 
frame. The channel that includes BCH is a single RF carrier frequency used for uplink 

10 and downlink. For high noise environments or for increased robustness, the BCH can 
hop frequencies according to a predetermined scheme or be repeated on several different 
frequencies. The repeating frame includes the downlink BCH for each base station, 
labeled BS1 etc. as shown in Table 3 below. The next frame includes the uplink 
Configuration Request CR, labeled CR1 etc. and downlink Configuration Message CM, 

15 labeled CM1 etc. Each frame also includes a number of reserved slots, shown as empty 
boxes below. These slots can be used for data traffic, if the broadcast channel is also 
used for traffic, for other control messages or reserved to reduce interference on other 
channels in the network. In one embodiment, the other traffic channels hop frequencies 
around and through the BCH. The frames are repeated for each respective base station 1 

20 to 128 to build a superframe as discussed in more detail below. After the last CM, 

CM 128, the superframe repeats and begins again with the next superframe and the BCH 
for base station 1 . 
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Uplink 


Downlink 


Superframe 
1 


Frame 1 








BS1 






Frame 2 


CR1 






CM1 






Frame 3 








BS2 






Frame 4 


CR2 






CM2 












Frame 255 








BS128 






Frame 256 


CR128 






CM128 






Superframe 
2 


Frame 1 








BS1 






Frame 2 


CR1 






CM1 













Table 3 



In another embodiment, the BCH is on its own channel and CR and CM are on a 
separate control channel. Alternately, one BCH can be provided on a constant frequency 
and a secondary BCH can be provided on another channel with hopping frequency. The 
hopping channel is described in the CM. 

Registration 

In one embodiment of the present invention, a session begins when a user 
terminal seeks to register with a base station. The user terminal does this without 
knowing its relative position to the best station in terms of direction and distance (or 
range) to the base station (BS). Therefore, when the user terminal (UT or remote terminal 
or subscriber station) requests registration using a configuration request (CR) burst, the 
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base station receives the transmission from the user terminal with a large timing 
uncertainty and unknown spatial parameters or weights when compared to traffic bursts 
that already utilize timing advance directives from the base station. 

Before registration, several user terminals may transmit CR bursts during the 
same time slot to register with the same base station. CR bursts addressed to different 
base stations may also be received. The base station resolves these requests using spatial 
processing and multiple antennas. By combining the multiple antenna measurements 
(beamforming) the base station minimizes the interference between signals and 
maximizes the signal-to-noise ratio (SNR) of each received burst. 

Beamforming can be performed in a variety of ways. In one embodiment, a 
training sequence is used with a least-squares cost function to determine the weights for a 
beamformer. This allows signals to be classified as desired or undesired. The training 
sequence-based approach uses an estimation of the timing and frequency offset of the CR 
burst. This estimation task can be performed by determining the least-squares cost for 
each timing and frequency hypothesis. This demands great computational resources. 
Using a special sequence design and a beamforming algorithm that does not require a 
search over the entire timing uncertainty range alleviates the computational load of the 
full search. In one embodiment, since CR bursts have a high timing inaccuracy, a 
periodic training sequence is dedicated for CR bursts alone. 

Configuration Request Burst Structure 

A configuration request (CR) burst is transmitted by a user terminal (user 
terminal) in order to initiate communications or registration with a base station (base 
station). It is transmitted after gathering information about the system by listening to 
(multiple) broadcast channel (BCH) bursts. The CR burst, is the first communication 
from a user terminal to the base station and therefore, the user terminal does not have any 
Docket No. 015685.P049 
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information about its range to its selected base station. Accordingly timing, range, and 
spatial processing weights among other things are all unknown to the base station. An 
example of a CR burst is shown in Figure 4. 

The configuration request burst is composed of several fields as shown in Figure 
4 which are listed in Table 4. The durations are described in terms of microseconds. In 
one embodiment, a symbol period is 2 microseconds. 



Duration 


Contents 


10 (isec 


ramp - up 


240 ^isec 


training symbols 


1 64 jasec 


information symbols 


10 (isec 


ramp - down 


1 06 (isec 


extra guard time 


15 (isec 


inter-burst guard time 



Table 4 Configuration Request (CR) Burst fields 



The training symbols are allocated 240 microseconds in order to allow the signal 
to be accurately received and demodulated before a user terminal has registered and has 
received any knowledge of the system. The training symbols are discussed in more detail 
below. 

The 82 information symbols are constructed from the configuration request 
message using, for example, forward error correction coding. In the present embodiment, 
the CR burst is modulated using n/2-BPSK modulation in order to decrease the peak-to- 
average ratio of the transmitted waveform. 
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The information symbols of the present CR burst are mapped out as shown in 
Table 5, below. Any of the items listed below can be deleted and transmitted later during 
the registration cycle or not at all based on the needs of the system. CR is scrambled by a 
function of BSCC ensuring that even if there is some interference from CRs sent to 
nearby base stations, the demodulation capture effect of the BSCC works out any 
collisions. In one embodiment, the scrambling is performed by taking the encoded bit 
sequence and exclusive OR'ing it with the output of a linear feedback shift register. 



Configuration Request Message 


Field 


# of Bits 


identity 


8 


utClass 


4 


txPwr 


5 


Total 


17 



Table 5 



identity is a set of unique random bits for each user terminal that differentiate 
simultaneous messages from multiple user terminals. Because of the randomness and 
large number of bits, it is unlikely that two user terminals will select the same identity 
code at the same time. 

utClass identifies user terminal capabilities (highest modulation class, frequency 
hopping capabilities, etc. ) This sequence identifies the type of user terminal that sent the 
CR. A palmtop digital assistant might have different capabilities than a desktop 
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computer with a fixed dedicated antenna. With utClass, the different capabilities can be 
distinguished. 

txPwr represents the power used by the user terminal to transmit the 
Configuration Request burst. For example, user terminal power = (2-txPwr - 30) dBm. 

In the present embodiment, the CR burst includes 106 microseconds of extra 
guard time. As mentioned above, the CR burst is transmitted by a user terminal without 
any knowledge of the time required for the burst to travel to the base station. In the 
present embodiment, the desired maximum range from a user terminal to a base station is 
15 km. The delay of the burst's reception at the base station includes the delay in 
receiving the BCH from the base station. A user terminal at range R km receives the 

o 

BCH burst(s) R/c seconds after they are transmitted by the base station (c=3 x 10 m/sec). 
When the user terminal decides to transmit a CR burst, its signal will be delayed by 
another R/c seconds, so that the delay of the CR burst at the base station relative to base 
station time will be 2R/c seconds. To accommodate delays created by a maximum range 
of 15 km, the guard period at the end of the CR burst time slot should be at least 100 
microseconds. In addition, in the present embodiment, the user terminal randomly delays 
its CR burst transmission from 0 to 9 symbols in order to reduce collisions with other 
user terminals that are sending CR bursts in the same slot and that may be at about the 
same range from the base station. Therefore, a CR burst from a distant user terminal may 
arrive at the base station as much as 1 18 microseconds late as compared to a user 
terminal that is very close to the base station and that does not delay its CR burst 
transmission. This accounts for the 106 microseconds of extra guard time that is allotted 
to the CR burst and the extra margin supplied by the 1 5 microseconds of inter-burst guard 
time. With the inter-burst guard time 3 microseconds of extra margin are allotted before 
the beginning of the next slot. 
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The present embodiment also includes an additional 15 microseconds of interburst 
guard time. This helps to confine the RF energy from the CR burst to the time-slot 
allocated for the CR burst in order not to create adjacent slot interference. The interburst 
guard time provides a total guard time of 121 microseconds to reduce the likelihood that 
5 a CR burst will extend into the next time slot. 

The CR burst is sent on the control carrier, as an example, exactly 2265 usee after 
receipt of a downlink BCH burst. In this way, an otherwise uninitialized user terminal can 
send CR without any knowledge of the frequency hopping sequence parameters. As 
discussed above, the CR burst is shorter than a standard uplink time-slot to allow for 
10 unknown time-of-flight from the user terminal to the base station and typically arrives 
late in the uplink time-slot receive window. 

Training Sequences 

1 5 In the present embodiment, known periodic training sequences are used for the 

training of base station beamformers in order to extract the desired CR bursts in the 
presence of interference. The training sequence includes several periodic repetitions of a 
core sequence followed by a marker sequence. Alternatively, the marker sequence may 
precede the repetition of the core sequence. By repeating the core sequence, the search 

20 range for a least squares (LS) beamformer can be reduced to a single repetition period of 
the training sequence. The period of the core sequence is a compromise between the 
duration of the search range and the possibility of finding sequences of the specified 
period with good autocorrelation and cross correlation properties. In the present 
embodiment, constraints on the periodicity, mean, autocorrelation and cross correlation 

25 are used to define the core sequences of the training sequence for CR burst acquisition. 
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Periodicity: x(k+P)=x(k), in which x(k) is the periodic part of a n/2-BPSK 
training sequence and P is the period of the sequence. In the present example, P=12 
symbols. 

Mean: the mean of the periodic sequence is constrained to be zero: 
p 

i Z x(k) = 0. 

* k=\ 

Autocorrelation: The out-of phase (symbol spaced) autocorrelation is bounded by 
1/3, where x is the phase offset of the out-of-phase sequence: 



I 



^ x * (k)x(k + r) 



i>(*)f 



<l/3, r*0 



Cross correlation: Symbol spaced cross correlation between two candidate 
periodic sequences x(k) and y(k) is bounded by 1/3: 

p 

k = \ 



*(k)y(k + r) 



At least two 7i/2-BPSK sequences conform to these constraints. More sequences 
may be possible if the constraints are relaxed. These two sequences are used as the core 
12 symbol sequence that is repeated to form the periodic training sequence as discussed 
below. They are: 

si = [1 j,l j,l,-j,-l j,-l,-j,-l,-j] 
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The bounds on auto and cross correlations help to make delayed versions of these 
sequences to appear partially uncorrelated to an LS beamformer which resolves them. 

The two periodic sequences can be shared among base stations. In one 
embodiment, depending on its color code, a base station will accept one of the two core 
training sequences but not both. The user terminal will be able to transmit the correct 
accepted sequence after learning the base station color code from the broadcast (BCH) 
burst. This enables a reuse pattern of 2. For reuse patterns greater than 2, some of the 
constraints on the sequences can be relaxed to yield a larger set of possible periodic 
sequences to distribute among the base stations. 

In an alternate embodiment, the user terminal selects which of the core sequences 
it will transmit. This can aid a base station in resolving two interfering CR bursts. The 
core sequence can be selected based on a serial number, product number, ID number or 
other stored number of the user terminal. For, example if the two sequences Si and S2 
above are available, the user terminal can read its serial number register. If the serial 
number is odd it selects sequence Si and if the serial number is even it selects S2. In 
another embodiment, the user terminal can generate a random number from 1 to 2. If 
there are more than two sequences a greater range of random numbers can be generated. 
The random number is then mapped into a list of sequences in order to select the periodic 
sequence that is associated with the generated random number. 

Transmitting CR Burst and Random Delay 

In one embodiment, the base station transmits broadcast channel (BCH) bursts 
continuously over a pre-established known channel. To begin a session with a desired 
base station, an unregistered and unrecognized user terminal in an unknown location 
listens to the BCH bursts. From the BCH, it determines the base station with which 
wants to establish communications and selects an appropriate training sequence for that 
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base station. It also selects a random number D between 0 and 9, and transmits its CR 
burst after delaying it for D symbols. The number for the delay can be selected as 
described above with respect to the periodic training sequence. The user terminal can 
select an internal register such as a Serial No. or Product ID No., that contains a number 
5 greater than 9 and use this as a pointer to select the number D. Alternatively, the user 
terminal can use a random number generator to generate a random number between 0 and 
9 and apply this number to select D. 

In one embodiment, the maximum value of random delay D is 9 symbols. 
Therefore, a user terminal can pick one of the 10 random delays, { 0,1, 2,. ..,9 } to send its 

10 CR burst to the base station. In addition, the maximum allowable roundtrip delay from 
the base station to the user terminal is assumed to be 50 symbols. Therefore, a CR burst 
can be delayed as much as 59 symbols for a user terminal at the maximum possible 
range. This delay is measured with respect to the CR burst of a user terminal that is at 
zero range and that did not delay its transmission. 

15 The base station receives and resolves the CR burst as described below. The 

roundtrip delay is determined to be 2(R/c) + D microseconds (since a symbol period is 2 
microseconds where R km is the range of the user terminal. The base station, in the 
Configuration Message (CM) burst, then instructs the user terminal to advance its time by 
Ta=(R/ c+ D/2) microseconds in its response to the CR burst. The base station is unaware 

20 of what, if any delay D the user terminal has applied to its CR burst. In response to this 
instruction, the user terminal, knowing its random delay D will subtract out the random 
delay from the propagation delay and advance its timing by x a -(D/2) microseconds. 

The random delay of D symbols is particularly valuable in a microcellular setting, 
or when user terminals cluster at about the same range from the base station. If the range 

25 is about the same for two different user terminals, the CR bursts from the two user 

terminals will look identical at the base station and may not be resolved. With random 
delays, the CR bursts will be spaced apart making them easier to distinguish. 
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CR Burst Resolution 

Figure 5 shows a training sequence of a CR burst for one embodiment of the 
invention. In Figure 5, x 510 is one of the two core sequences {si, S2}. x is the vector 
5 consisting of 12 symbols X] to X12. In this application, boldface text is used to represent 
vectors and normal text is used to represent scalars. The complete CR Burst training 
sequence is constructed from the core periodic sequence x. Nine repetitions of x; 
P1,P2,... 5 P9 are followed by a marker sequence M 505. Figure 5 also can be viewed 
from the perspective of base station timing. In one embodiment, Figure 5 shows 

10 windows that the base station uses to resolve received training sequences in a CR burst. 
In this context, the sequences PI to M would line up with the windows as shown only if 
the sequences were received from a user terminal that was in the same location as the 
base station, or in other words, if there were no propagation delay between the base 
station and the user terminal. In this paradigm, a round trip delay margin window 520 

15 contains the sequences PI through P5 for a total of 60 symbols or 120 microseconds. 
This accommodates a combined delay of 50 symbols due to the round trip propagation 
delay between the user terminal and the base station and 1 0 symbols of user terminal 
selected random delay. A beamformer analysis window 530 contains core sequences P6 
through P9 and so is 48 symbols wide. This window will capture 4 periods of x even if 

20 the CR burst is delayed by up to 60 symbols. 

Using the observations that fall into the beamformer analysis window, a base 
station can estimate the delay encountered by the CR burst modulo- 12 symbols since 
each period PI through P9 is identical and 12 symbols long. In other words, the 
measurement in the beamformer analysis window will be the same for the delays x and x 

25 +12, where 0 < x < 12. The marker sequence M at the end of the burst can be used to 
resolve this modulo- 12 ambiguity. 
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In order to demodulate the payload part of the burst, the absolute timing 
information instead of the relative timing or modulo- 12 timing information is desired. As 
mentioned above, the marker sequence (M) at the end of the training sequence of Figure 
5 can be used for this task. The marker sequence M is set to be M = - x so that one period 
of the core training sequence x augmented with M will be orthogonal to two concatenated 
periods of x, i.e., [x, M] [x,x] H = 0 in which x = [xi,..., X12], the core sequence, forms a 
period of the training sequence. 

Since, in the present embodiment, the maximum allowed delay is 60 symbols, a 
roundtrip delay margin window with a duration of 5 periods of the core sequence is 
provided. This leaves a duration of 4 periods of the core sequence for the beamformer 
analysis window. Absolute timing is determined by determining which of the four 
possible periods of the core sequence is the last one received before the marker sequence. 
The beamformer analysis window may contain portions of five different core sequences, 
at least three of which are full periods of the core sequence. The first symbol of the 
found sequence is either 12, 24, 36 or 48 symbols from the first symbol of the marker 
sequence. Therefore the absolute timing can be determined by testing only five 
hypotheses. In other words, given x, (0 < x < 12), as the modulo- 12 delay estimate from 
the spatial processing search, the possible candidates for the absolute delay are: {x, x +12 
,x+24,x + 36,x+48 }. 

As can be seen from Figure 5, for any delay less than 60 symbols, the beamformer 
analysis window at the base station will always contain at least 4 periods of the periodic 
sequence x. As a result, the first symbol, x\, of the first repetition of the core sequence 
within the beamforming analysis window can be at most 12 symbols away from the 
border of the round trip delay margin window and the beamformer analysis window. 

This property of the repeating periodic training sequence can be taken advantage 
of when searching for the first symbol, xi, of the core sequence x. Figure 6 shows a flow 
diagram for resolving the CR burst according to one embodiment of the invention. The 
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CR burst is received at the antenna array 602 and the measurements from the multiple 
antennas for the period that covers the beam former analysis window are stored 604. 
Then the beamforming weights that will yield the minimum least squares (LS) error for 
each timing hypothesis in the search window of 12 symbols is determined. The least 
5 squares error results are used to find a coarse relative timing hypothesis 606. Finding this 
timing includes performing a coarse search for the symbol timing modulo- 12 and 
performing a follow-up search 608. 

A weight vector is then determined 610 using the relative timing. This weight 
vector is applied to the stored measurements to convert the measurements from each 
10 antenna channel into a single channel 612. Then, this single channel is analyzed to 
g determine the fine timing 614 and the frequency offset 616 of the desired signal . The 

fine timing and frequency estimates are used to determine a new more accurate weight 
vector 618 based on the data in the beamforming analysis window. This improves the 
quality of the weight vector and the new weight vector is applied to the stored 
1 5 measurements to convert the measurements from each antenna channel into a new single . 
channel 620, further improving the quality of the received signal. The processing so far 
cannot perfectly determine the timing information, because it can only measure time 
delays modulo-12 symbols, the duration of the repeating core sequence. The timing 
estimate is in the range {0... 12}, i.e., delays of x and x +12 symbols will ideally produce 
20 the same timing estimate at this point. Using the new weight vector, absolute timing is 
determined by finding the marker sequence 622. This resolves the modulo-12 timing 
ambiguity. 

When the aBSoJute timing is known, the entire training sequence can be identified. 
Accordingly a new weight v&qtor calculation can be made using the entire training 
25 sequence including the marker seqbsrure 624. This redetermined weight vector is still 
more accurate because of the larger numb^isof samples that can be used. It is applied to 
the stored measurements to again convert the me^jjrements from each antenna channel 
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into a single channeh^26. With the new single channel, the CR burst is demodulated and 
read 628. The redeterminetk^veight vector is also used to send a CM (configuration 
message) burst or other traffic bacKto the user terminal 630. 

Thisfcocessing structure yields accurate beamformer weights that can be used to 
receive the desirecKignal without determining perfect timing for the signal. The 
beamformer weights caii^lso be used to transmit a burst with appropriate spatial 
directionality and interferenc^nd noise suppression back towards the user terminal in 
the appropriate corresponding dowhHnk slot. 

Coarse Search for Relative Timing 

After the CR burst has been received in each of the antennas of the array 602 and 
the signals have been stored 604. The system can proceed to determine the relative 
timing based on the stored measurements 606. In one embodiment, this process is done 
using a covariance matrix computation and a Cholesky decomposition. 

For measurements that lie in the beamforming analysis window assume that 
measurements are collected at 1 .5 times the symbol rate. In this example, the 
beamforming analysis window is of duration 48 symbols or 4 repetitions or periods of the 
core sequence. Next, assume that the variable coreSnapPoint points to the first snapshot 
that is in the beamformer analysis window, and r(t) denotes the snapshot (measurement 
vector) collected at time t seconds after the start of the CR burst (relative to base station 
time). The covariance matrix R can be estimated as follows, where Tbaw is the time of 
the beamforming analysis window: 

R= £r(t)r H (t), 

Tbaw ~ {'| coreSnapPoint < / < coreSnapPoint + 48 } 

te T BAW 

After computing the covariance matrix R, a matrix L can be found such that 
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R = L L H 



L is termed as the Cholesky factor of R. 

In one embodiment, the location of the first symbol of the core training sequence 
over the first 12 symbols of the beamforming search window is searched for using a 
least-squares processor with a search step of 2/3 symbols. This results in 18 possible 
locations of the first symbol, i.e. 18 hypotheses and, accordingly, 18 least-squares (LS) 
error calculations. The vector coarseSearchGrid contains the delay hypothesis values 
for each 2/3 symbol increment (The unit used here is a symbol period): 

coarseSearchGrid = [0,2/3,4/3,2, ...,1 1 1/3] which contains the delay values for 
the 1 8 hypotheses. 

For each delay hypothesis k, (1 < k <18), the following computations are 
performed: 

Calculate cross-correlation vector pk for each hypothesis: 
Pk= £ r(/)rf *(f+r 4 ) 

In which r(t) is the received signal samples and d(t) is the samples of the desired 
signal. The desired signal samples are assumed to be available in a filter bank over- 
sampled at 24 times the symbol rate. The value ik is the assumed delay for the k-th 
hypothesis, i.e., ik = coarseSearchGrid (k). 
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Back Substitution: The resulting vector pk is applied to solve for an interim vector 
Xk for each hypothesis k, where L is the Cholesky factor as mentioned above: 

Lx k = p k 

5 

Then the least-squares (LS) fit for each hypothesis is computed: 
f k = f(h)= Z \d(t+r k )\ 2 -Xk H X k 

Due to the periodicity of the desired signal within the beamforming analysis 

window, the cross correlation vector computation can be modified to: 
10 P* = X l>(0+ r(f + 12)+ r(f + 24)+ r(f + 36)]d * (f + r k ) 

in which Ti is a period of the desired signal, i.e., 
Tj = { t |coreSnapPoint < t < coreSnapPoint +12} 

Therefore, the above calculation of the cross correlation vector Pk costs one-fourth 
15 as many multiplies as it would if the core sequence were not a repeated periodic sequence 
within the training sequence. Finally, the hypothesis with the minimal LS error f k is 
picked. An illustration is presented with respect to Figure 7. 

Figure 7 shows the LS cost function on the vertical axis 702 as a function of the 
time delay x for the 1 8 hypotheses on the horizontal axis 704. The horizontal axis is 
20 marked in units of 2/3 of a symbol time in accord with the sampling rate of 1 .5 times the 
symbol rate. The selected best coarse search time delay estimate 706 is shown in Figure 
7 at t = 4/3 symbols. Note that the delay is measured with respect to the nominal time or 
modulo- 12 time at which the CR burst is expected to arrive. This nominal time is named 
the coreSnapPoint. This simplified computation of the cross correlation vectors pk based 
25 on the periodicity of the training sequence within the beamforming analysis window 
yields substantially the same accuracy as a regular cross correlation computation. Any 
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mismatch of the training sequence in this simplified approach makes only a minimal 
impact on the performance of the beamformers due to the periodicity of the training 
sequence. 

Follow-up Search 

Having made a coarse determination of the relative timing 606 a follow-up timing 
search 608 can be performed. This search has one more LS error calculation. The result 
of this calculation can be compared to the best LS error from the coarse search. 

The LS error vector f k from the coarse search is treated as a circular buffer and 
concentration is on the midpoint of the range between the best location from the coarse 
search (that yielded the minimum LS error) and its neighbor that yielded the next best LS 
error. 

In the coarse search, the following 20 dimensional vectors were formed : 

updatedLSValues = [fi 8s fi,f2,f3,f4, , fn,fi8, fi] 
and 

updatedCoarseSearchGrid = [-2/3,0,2/3,4/3,3 1/3, 4, ...,11 1/3, 12] 

Assume that the hypothesis of f3 (x = 4/3) (of the current list) was selected by the 
coarse search procedure (as shown in Figure 5) because f3 < f2 < f4. In the follow-up 
search procedure, the LS error is computed for x = 1, since this forms the midpoint of the 
range between the delays corresponding to fi (x = 2/3) and f3 (x = 4/3). This is shown 
graphically in Figure 8. 

Figure 8 shows the LS cost function 802 as a function of the 1 8 hypotheses 804 in 
which the time domain is again sampled at a rate of 1 .5 times the symbol rate. The 
selected best case coarse search time delay estimate 806 is x = 4/3 symbols. After the 
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additional hypothesis test, the hypothesis (among the 19 of them) which corresponds to 
the minimum least-squares error is again picked. In this example, x 0 = 1 is the delay 
value 808 that accomplishes this. 

Beamforming Weights 

After the follow-up search, determining the relative timing delay x to within one 
third of a symbol period 608, the weight vector w 0 is recomputed to extract the desired 
signal as follows: 

Using the cross-correlation vector p 0 for the best hypothesis: 
p 0 = £r(0</*(/ + O 

And the corresponding interim vector x 0 : 
Lx 0 = p 0 

The weights can be computed by solving the following back substitution: 
L H w 0 = p 0 

With the resulting LS error which is already determined to be f(x 0 ): 

f(Xo)= X |d(t + X 0 )| 2 -X 0 H X 0 

The beamforming operation results in a single channel measurement that shall be 
named g(t) and defined as follows: 
g(t) = w 0 H r(t) 

The actual computations can be simplified because the time range of the 
beamforming operation does not need to be the whole burst. Note that the frequency 
estimation (search) has not yet been performed. This reduces the search to a single 
dimension, namely time-delay. The frequency search can be avoided because of the 
duration of the beamforming analysis window, the maximum frequency deviation 
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expected in the received signal and the target SNR level at the output of the beamformer. 
When the beamformer analysis window is short enough to make the frequency offset 
negligible, the desired signal mismatch is not significant enough to degrade performance. 
Furthermore, the inability to improve the SNR even by using an MSE beamformer (that 
5 knows the perfect covariance matrix and steering vector) to high levels (above 5 dB) 
makes the actual beamformer blind to imperfections in the desired signal (such as 
frequency offset). The presence of frequency offset becomes serious at high target SNRs 
and ignoring large frequency offsets causes the SNR performance to tend to saturate at 
around 5 dB output SNR. However, this is still not an issue if the CR burst payload can 
10 withstand SNRs as low as 0 dB. Furthermore, the single channel estimation of frequency 
and fine timing that follows will alleviate these problems using fewer computations. 

Fine Timing Estimator 

1 5 The coarse search and the follow-up search provide a timing delay with an 

effective rate of three times the symbol rate. The timing resolution can be improved to 
better demodulate the payload. Severe intersymbol interference can be expected at the 
worst case timing events (1/6 symbol period of timing error). To estimate fine timing 
614, the Oerder-Meyr blind timing estimator can be used. This estimator is blind in the 

20 sense that it does not assume any knowledge of the transmitted symbols. Therefore it is 
insensitive to frequency offset errors. The Oerder-Meyr estimator can be used as follows: 

First, interpolate g(t) to 3 times the symbol rate (from 1 .5 times the symbol rate) 
using an interpolation filter, and call the output gj(t). The duration of the measurements 
process amounts to 48+12+82=144 symbol periods, this is the guaranteed maximum 

25 number of symbols on the right hand side of the coarse timing estimate in the 

beamforming analysis window as shown in Figure 5. Alternatively, the symbols can be 
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sampled at a three times rate originally. In such a case, only every other sample would be 
used in the coarse search and follow-up search. 

Next, the interpolated signal is passed through a memoryless nonlinearity to 

obtain 



The absolute value nonlinearity generates tones at the multiples of the symbol 
rate. The nonlinearity is approximated by fitting a polygon. For example, if z is a 
complex number with representation z = z r + jzi, then \z\ « max(|z r |,|zi|) + 0.34 x min 
(|z r |,|zj|). The zero-crossing of the complex sine wave at the symbol rate is determined by 
10 using a 3 point DFT and this location is termed as a symbol transmission instant. Finally, 
the symbol transmission instant is found that is closest to the position identified by the 



follow-up search as the location for the first symbol of the periodic training sequence. 

The closed-form estimate obtained using the Oerder-Meyr timing estimator is 
termed as Tn n e. If this estimate is not in the range (0,12), it is mapped into this range by 
1 5 adding or subtracting 12 to it. 



20 window together with the desired signal to compute the cross correlation function and 

determine the coarse frequency estimate 616 as the peak of the cross correlation function. 
For simplicity, this analysis is done at the 1.5 times over sampling rate (instead of the 3 
times rate used for the fine timing search), however many other sampling rates can be 
used. Specifically, for the m-th frequency hypothesis, the function freqMatch(m) is 

25 computed : 



5 



gn(0= |gi(t)|. 



Coarse Frequency Estimator 



After fine timing is acquired 614, the measurements are used in the analysis 



freqMatch(m) = £ g{t)d * (/ + r fme ) exp( -jlnfj) 
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The frequency candidate that maximizes freqMatch(m) is selected as the 
frequency estimate, f CO arse- Frequency candidates (f m 's) are in the set: 

{-2800,-2000,-1200,-400,400, 1200, 2000,2800} 

5 

Alternatively, the frequency uncertainty can be regarded as much less and fewer 
hypothesis tests can be performed. 

More Accurate Beamforming 

10 

Using the frequency and timing estimates Tf ine? fcoarse, an improved weight vector, 
Wfjnai can be determined 618 and applied to the found core sequences to obtain an 
estimate of the desired burst S(t): 

15 s(t) = w fina , H r(t) 

The more accurate weight vector is obtained by employing d(t+Tf, ne ) exp (j27iffj ne t) 
as the desired signal in the LS error computations. However, if the corresponding LS 
error is larger than the follow up search based LS error f 0 , the fine timing and frequency 
20 estimates are rejected and the estimates are set to Tfi ne = x 0 and f CO arse = 0, which results in 

Wfi na | = W 0 . 

To further improve the SNR results at high post-copy SNR levels, the search can 
be performed for a fine frequency estimate around the coarse frequency estimate at 
locations f CO arse + [-800,0,800] using an LS fit to multichannel measurements for t e 
25 T B aw- Such a search can be included to improve the frequency estimation performance 
for the CR burst demodulator. 
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Absolute Timing Recovery 

Prior to determining the absolute timing, as opposed to relative timing, the more 
accurate weight vector is applied in order to convert the measurements stored from each 
5 of the antennas into a single channel measurement 620. This conversion is still done 
using the data in the beamforming analysis window because only this data is certain to 
consist only of periods of the core sequence. An appropriate procedure 622 for 
determining the absolute delay can be implemented as follows. Assume that the modulo- 
12 time-delay estimate is x, and that g(t) denotes the output of a beamformer (as 
10 described above g(t) is a single channel measurement) that extracts the desired CR burst. 
First, the candidate starting locations for the concatenated sequence [P9,M] (see Figure 5) 
for each of the five delay hypotheses is determined. For example, if the absolute delay 
equals x, the sequence [P9,M] should start x + 36 symbols away from the left border of 
the beamformer analysis window. For the hypothesis in which the absolute delay is 
15 assumed to be equal to x + 48, the sequence [P9, M] should start x 84 symbols away from 
the left border of the analysis window. 

More specifically, assume that Tk denotes the begin time for the k-th hypothesis, 
i.e., Tk = x+36+(12xk) 5 for k=0,l ,2,3,4. Let \u denote the samples of g(t) that correspond 
to the interval which starts at Tk for a duration of 24 symbols (two periods). For each 
20 hypothesis, compute: 

h(k) = |v k [x,-x] H | 2 - 

Pick the hypothesis that yields the maximum h(k). 

Assume that hypothesis ko is selected. Then, the absolute timing estimate is 
determined as x + (koX 12). 
25 Due to the orthogonal nature of [x, m] and [x, x], there is a five-ary pulse-position 

modulation (PPM) demodulation problem in which the SNR changes depending on the 
true position of the desired vector [x, m]. For example, if the actual delay is less than a 
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period, then the payload portion of the desired CR burst will appear as noise for the 
hypothesis tests 1 ,2,3,4. On the other hand, if the actual delay is, for example 20 
symbols, then, the CR burst payload symbols will not contribute to the test of a 
hypothesis in which k=0 as noise because of the orthogonality of [x 5 m] and [x, x], but 
5 they will contribute to the noise component in h(k) for k=2,3,4. With this delay condition, 
the beamforming analysis window receives the following 48 symbols: 

X 5 ,X6 ? X7,X 85 X9,Xio 5 Xii,Xi2,X],X2,X3,X 4 , ... X 5 ,X 6 ,X7 3 X 8 ,X 9 ,Xio,Xi i,Xi2,Xi,X2,X 3 ,X4 9 

In this stream, the first and last 1 2 symbols are written out. The spatial processing 
search will discover the first occurrence of xi at a delay of 8 symbols. The next step is to 
10 determine whether the absolute delay is 8,20,32,44,56 symbols. In order to determine 
this, the measurements following the beamforming analysis window are focused upon. 
With a delay of 20 symbols, the next 40 symbols will be: 

Xs^Xy^X^XiOjXi^X^^i^^^Xs, X 6 ,X7,X 8 ,X 9 ,Xio,Xl |,X| 2 ,-X] ,-X2,-*3,-X4, 

-x 5 ,-x 6 , -x 7 ,-X8,-X9rX i or* 1 1 >-x 1 2,di ,d2,d3,d4 5 d 5 ,d6,d7,d8 
1 5 in which dj indicates the i-th payload symbol, which is unknown. As a result, vo = 

[x, x] and vi = [x, -x], which results in h(0)=0 and h(l) = 24 2 . Accordingly, the desired 
signal does not contribute to the decision variable h(0) when the absolute delay is larger 
than 12. On the other hand, v 2 = [-xi,-x 2 , ... ,-Xi2, di,d2,...di2] is not necessarily equal to 
zero. The first half of \2 equals the marker sequence, which yields the result: 
20 h(2) = |-12-[d,,d 25 ...d 12 ]x H | 2 . 

This differs significantly from: 

h(3)=|[ di,d 2 ,...d 12 ] x H - [d, 3 ,d,4,..., d 24 ] x H | 2 . 

h(2) is expected to be higher than h(3) on average because of the coherence term 
25 (-12) inside the magnitude operator and its presence creates a significant competition to 
the correct hypothesis decision value (h(0)) in the presence of measurement noise. 
Therefore, absolute delay errors that are equal to +12 should dominate other values of 
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delay estimation errors. One way to eliminate the self-noise problem is to put the marker 
sequence in the beginning of the burst and search for the concatenated sequence [-x, x]. 
In such a design, self noise will exist only due to the presence of the ramp-up symbols, 
but this will affect the noise component in only one of the h(k)'s. Other modifications to 
5 the ordering of the sequences are also possible to suit particular applications. 

As discussed above, the remaining task is to find the marker sequence. This can 
be done based on a cross correlation operation with the 1 .5 times over-sampled version of 
the sequence [x,-x] with Vk corresponding to all of the five hypotheses locations 
indicating period slips of 0,1,2,3,4. 
1 0 Using the fine timing estimate Tf me? the time of the final beamformer output 

sample (that arrives at t c ) which is closest to coreSnapPoint + Tfj ne is determined. Let x e 
denote the time difference between the chosen measurement at t c and the location of 
actual timing location: 

15 x e = t c - coreSnapPoint -ifme 

Remember that the beginning of the k-th hypothesis (k=0,l ,2,3,4) will be at times 
T k = coreSnapPoint + Tf ine . Therefore, the vector v k to be tested for the cross correlation 
for the k-th hypothesis will take the form: 



v k = [S(tc + 12 x (3 + k)) , S(tc + 12 x (3 + k)+ 2/3)], 



s(tc + 12 x (3 + k)+ 2 x 2/3 ), ... , s(t c + 12 x (3 + k)+ 35 x 2/3 ) 



which is a vector with 36 entries. 



The v k 's are cross correlated with the desired signal vector d, which is formed 



from the knowledge of the training waveform and the timing difference i e : 



Docket No. 015685.P049 

Express Mail No. EL609099423US 



38 





d = [Sd(coreSnapPoint+ x e + 12 x 3), Sd(coreSnapPoint+x e +12 x3+2/3), 



Sa (coreSnapPointH- x e + 12 x 3+ 35 x 2/3) ] 



in which Sd(t) is the training waveform for the desired signal. 



The magnitude squared cross correlations for each hypothesis is computed: 



h(k) = |v k d H | 2 



Let ko denote the index of the largest h(k). Then, the absolute timing is 
determined as: 

Tabs = Xfine + (koXl2) 

The values are in terms of symbol durations in the above equation with the period 
of the training sequence being 12 symbols. 

Alternatively, symbol spaced sampling can be used in the absolute timing 
recovery. This can eliminate the mismatches created in the 1 .5 times measurements that 
result in loss of orthogonality. With symbol spaced sampling, this operation will require 
5 cross correlations of length 24 and comparing the squared-magnitudes of the cross 
correlations. 

After absolute timing recovery 622, the whole training sequence can be used to 
improve the timing and frequency estimates and these estimates can be used to perform 
one more stage of beamforming over the whole training sequence duration. This can 
increase the number of symbols used in spatial processing to 120 symbols from 48 and 
improve the performance at a cost of computational complexity The beamforming would 
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involve first determining a still more accurate weight vector using the full training 
sequence 626 and then applying this weight vector to convert the measurements stored 
for each antenna channel into a single channel measurement 628. In doing so a new 
covariance matrix should be computed and therefore another Cholesky factorization 
5 performed. The conversion can be applied to the entire CR burst so that the CR burst can 
be demodulated as a single channel 628. Alternatively, the absolute timing can be 
applied to determine the beginning and end of the CR burst and the existing weight 
vector and frequency estimates can be used to demodulate the signal. The choice will 
depend on the demands of the particular system and the condition of the channel. Having 
10 determined timing and a weight vector, this weight vector can be used to send traffic back 
to the remote from which the CR burst was received 630. 

General Matters 



15 The present invention can eliminate the need to acquire the absolute timing 

information in order to perform beamforming. It exploits the periodicity of target 
training sequences in order to estimate the time delay encountered by the desired signal 
modulo-the period of the training sequence which is less than the total time delay 
ambiguity range. The reception of uplink messages at the base station which are not 

20 perfectly synchronized is enhanced. 

Training based interference cancellation methods can be used to enhance the 
signal quality for such bursts. However, the training sequence should be aligned with the 
incoming signal, and therefore the resulting search process involves testing every 
possible delay with a predetermined resolution. As the delay ambiguity grows, this 

25 hypothesis testing process becomes increasingly more expensive with computational 
resources. 
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Using a periodic training sequence eases the efforts of hypothesis testing followed 
by a timing acquisition sequence to determine the perfect timing information. The 
periodic sequence decreases the timing ambiguity to within a period of the training 
sequence, which can be selected small so that resources are available to handle the 
5 hypothesis testing task. For example, if the expected delay is 50 symbols, and if the 
period of the training sequence is 10 symbols, with a resolution of a symbol, 10 
hypotheses can be tested instead of 50, provided that the periodic training sequence can 
be observed in the search window, regardless of the delay. In addition, the search 
operation, which involves cross correlations can be simplified by the number of periods 

10 in the search window, if the search window contains 5 periods of the training sequence 
(period with 10), the cross correlations can be simplified by five-fold. 

Once the timing hypothesis is verified, only the number of periods of the 
repeating core sequence that has been skipped must be determined. To perform this more 
accurately, the multiple channel measurements (if available) can be converted to a single 

15 channel measurement using the search procedure results outlined above. After that a 

timing acquisition sequence that is at the beginning of the burst can be found, and to find 
the number of periods that were missed. So in the previous example, only periods 
{0,1,2,3,4} may have been missed in the beamforming analysis window. Accordingly, 5 
searches can be performed and these are performed on single channel measurements 

20 rather than multichannel measurements. 

The invention eliminates the high cost of search on possibly multichannel data for 
all possible delays that may be encountered by the user terminal as a function of range. 
The burst structure and periodic training sequence makes it easier to search for the arrival 
time of the signal and the periodicity of the training further reduces the cost of 

25 computations. 

In the description above, for the purposes of explanation, numerous specific 
details are set forth in order to provide a thorough understanding of the present invention. 
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It will be apparent, however, to one skilled in the art that the present invention may be 
practiced without some of these specific details. In other instances, well-known structures 
and devices are shown in block diagram form. 

The present invention includes various steps. The steps of the present invention 
5 may be performed by hardware components, such as those shown in Figures 1 and 2, or 
may be embodied in machine-executable instructions, which may be used to cause a 
general-purpose or special-purpose processor or logic circuits programmed with the 
instructions to perform the steps. Alternatively, the steps may be performed by a 
combination of hardware and software. The steps have been described as being 

10 performed by either the base station or the user terminal. However, any steps described 
as being performed by the base station may be performed by the user terminal and vice 
versa. The invention is equally applicable to systems in which terminals communicate 
with each other without either one being designated as a base station, a user terminal, a 
remote terminal or a subscriber station. 

1 5 The present invention may be provided as a computer program product which 

may include a machine-readable medium having stored thereon instructions which may 
be used to program a computer (or other electronic devices) to perform a process 
according to the present invention. The machine-readable medium may include, but is 
not limited to, floppy diskettes, optical disks, CD-ROMs, and magneto-optical disks, 

20 ROMs, RAMs, EPROMs, EEPROMs, magnet or optical cards, flash memory, or other 
type of media / machine-readable medium suitable for storing electronic instructions. 
Moreover, the present invention may also be downloaded as a computer program product, 
wherein the program may be transferred from a remote computer to a requesting 
computer by way of data signals embodied in a carrier wave or other propagation 

25 medium via a communication link (e.g., a modem or network connection). 

Importantly, while the present invention has been described in the context of a 
wireless internet data system for portable handsets, it can be applied to a wide variety of 
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different wireless systems in which data are exchanged. Such systems include voice, 
video, music, broadcast and other types of data systems without external connections. 
The present invention can be applied to fixed remote terminals as well as to low and high 
mobility terminals. Many of the methods are described in their most basic form but steps 
can be added to or deleted from any of the methods and information can be added or 
subtracted from any of the described messages without departing from the basic scope of 
the present invention. It will be apparent to those skilled in the art that many further 
modifications and adaptations can be made. The particular embodiments are not 
provided to limit the invention but to illustrate it. The scope of the present invention is 
not to be determined by the specific examples provided above but only by the claims 
below. 
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